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ABSTRACT 



Recently, high-dispersion spectroscopy has demonstrated conclusively that four of the five globular clusters (GCs) in the Fornax 
dwarf spheroidal galaxy are very metal-poor with [Fe/H] < -2. The remaining cluster, Fornax 4, has [Fe/H] = -1.4. This is in stark 
contrast to the field star metallicity distribution which shows a broad peak around [Fe/H] » -1 with only a few percent of the stars 
having [Fe/H] < -2. If we only consider stars and clusters with [Fe/H] < -2 we thus find an extremely high GC specific frequency, 
5,v ~ 400, implying by far the highest ratio of GCs to field stars known anywhere. We estimate that about 1/5-1/4 of all stars in the 
Fornax dSph with [Fe/H] < -2 belong to the four most metal-poor GCs. These GCs could, therefore, at most have been a factor of 
4-5 more massive initially. Yet, the Fornax GCs appear to share the same anomalous chemical abundance patterns known from Milky 
Way GCs, commonly attributed to the presence of multiple stellar generations within the clusters. The extreme ratio of metal-poor 
GC- versus field stars in the Fornax dSph is difficult to reconcile with scenarios for self-enrichment and early evolution of GCs in 
which a large fraction (90%-95%) of the first-generation stars have been lost. It also suggests that the GCs may not have formed as 
part of a larger population of now disrupted clusters with an initial power-law mass distribution. The Fornax dSph may be a rosetta 
stone for constraining theories of the formation, self-enrichment and early dynamical evolution of star clusters. 
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1. Introduction 

The Fornax dwarf spheroidal galaxy (dSph) is well-known 
for its high globular cluster (GC) specific frequency, i.e., the 
number of GCs normalized to a host galaxy My = -15 
(Harris & van den Bergh 1981). Assuming My = -13.2 (iMateol 
[l998), the five GCs (Hodge 1961) correspond to a specific fre- 
quency of Sn = 26. Specific frequencies rivalling that of Fornax 
have only been found in other similarly faint dwarf galaxies 
dPeng et al.l2008tlGeorgiev et al.l2010i) . Spiral galaxies typically 
have Sff ~ while normal elliptical galaxi es have ^ 3-5, 
slightly higher in clusters than in the field (lHarrislll99lh. Even 
the rnost cluster-rich cD galax ies generally have Sn ^ 15 dHarrisI 

1199 It iBrodie & Stradel^l2006^ ■ These differences constitute the 
classical "specific frequency problem". 

A second specific frequency problem becomes apparent 
when considering the numbers of GCs associated with differ- 
ent stellar populations within galaxies. Early-type galaxies are 
generally redder (more metal-rich) than the ave r age of their GC 
system s dForte et all 119811: iLarsen et al.1 120011; iForbes & Foitg 

1200 lb . and direct comparisons of metallicity distributions for 
stars and globular clusters have confirmed that the specific fre- 
quency t ends to increase with decreasing meta llicity within 
galaxies dHarris & Harrisll200l iHarris et al.ll2007h . This is also 
evident in our own Galaxy: of the roughly 150 GCs known in the 
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Milky Way, about 2/3 are associated with the halo dZinnlll985h . 
even though the halo only contains 1 %-2% o f the stellar mass 
dSuntzeff et anil99ll:lDehnen & Binnevlll998h . 

It has, in fact, been suggested that a significant fraction of 
the Milky Way halo stars originate from disrupted GCs. If the 
Galactic GC population formed with a power-law mass function 
similar to that obse rved in young cluster systems at the present 
epoch ('Larsenl l2009l) . then the stars lost from dissolving clusters 
over a Hubble time might be suffi cient to account for e ssentially 
the entire stellar halo ('e.g. lKruiissen & Portegies Zwartll2009h . 

The above estimate only involves standard dynamical evo- 
lution of the clusters. However, a similar conclusion can be 
reached via a different line of reasoning, starting from the in- 
creasing body o f evidence that GCs host dual or multiple stel- 
lar generations dGratton et al.l l2012h . In the following we will 
refer to the "first" and "second" generation of stars, keep- 
ing in mind that the actual star formation histories may be 
more complex. The second generation is characterized, among 
other things, by anomalous abundances of several light ele- 
ments that imply p-capture nucleosynthesis at high temper- 
atures. The cuiTently favoured sites are massive AGB stars 
(D'Ercole et al. 2008), fast-rotating massive main sequence 
stars ([Decressin et al.l 12007), or massive interacting binaries 
dde Mink et al. l2009h . A fundamental challenge faced by most 
scenarios for the origin of multiple stellar generations in GCs is 
the "mass budget" problem: at the present epoch, the second gen- 
eration of stars has a mass similar to, or even greater than the first 
generation. However, for a standard IMF the ejecta produced 
by the first-generation stars fall short by large factors compared 
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Table 1. Metallicities and absolute magnitudes for the Fornax 
GCs. 





[Fe/H] 


Mv 


Fornax 1 


-2.5 ±0.1' 


-52' 


Fornax 2 


-2.1 ±0.1' 


-13' 


Fornax 3 


-2.3 ±0.1^ 


-9,2' 


Fornax 4 


-1.4 ±0.1- 


-12' 


Fornax 5 


-2.1 ±0.1^ 


-lA' 



Re ferences. (1) iLetarte et alj j2006h : (2) iLarsen et al.l j2012D : 
(3) IWebbin]d (119851) 



to the mass required to form the observed second generation, 
even if a 100% star formation efficiency is assumed. This has 
led to the suggestion that a large fraction of the first generation, 
up to 90-9 5%, was los t soon after the formation of the second 
generation ( DlErcole et al. 2008; Schaerer & Charbonnel 20lll 
lBekldll20Tl ). Accounting for subsequent dynamical evolution, 
one again finds that a significant fraction of the Milky Way stel- 
lar halo might come from GCs, even without considering that a 
significantly larger population of lo wer-m ass clusters may have 
been present initially ( Gratton et al.ll2012l) . 

In this letter we discuss the globular cluster system of the 
Fornax dSph in the context of the second specific frequency 
problem and self-enrichment scenarios. While these two issues 
may seem unrelated, the Fornax dSph turns out to be such an ex- 
treme case of the second Sf^ problem that it puts tight constraints 
on the amount of mass that could have been lost from its GCs. 



2. Chemical composition of the Fornax GCs 

We have recently carried out a detailed analysis of the chem- 
ical composition of the three clusters Fornax 3, 4 and 5 from 
integrated-light high-dispersion spectroscopy with UVES on the 
ESO Very Large Telescope. O ur analysis techniq ues are de- 
scribed in a companion paper (iLarsen et al.ll2012l submitted), 
along with a more detailed discussion of the abundance measure- 
ments and references to previous determinations of the metal- 
licities. We measured iron abundances of [Fe/H] = -2.3 ± 0.1 
(Fornax 3), - 1 .4 + 0. 1 (Fornax 4) and -2. 1 ±0.1 (Fornax 5). Th e 
clusters Fornax 1, 2 and 3 were observed bv lLetarte et al.l(l2006h . 
who measured abundances for three individual RGB stars in 
each cluster; they obtained [Fe/H] = -2.5 ± 0.1 (Fornax 1), 
-2.1 + 0.1 (Fornax 2) and -2.4 ± 0.1 (Fornax 3). We note that 
there is excellent agreement between the two independent mea- 
surements for Fornax 3. The high-dispersion measurements of 
the Fe abundances are summarized in Table [1] It is clear that 
Fornax 4 is significantly more metal-rich than the other clus- 
ters, all of which have [Fe/H] < -2. We also list t he GC My 
magnitudes, based on the magnitudes in IWebbinkI (fT98l and 
using the same apparent distance modulus of m - M - 20.80 
as for the Forna x dSph itself for consistency dMateol Il998t 
iPietrzvriski e"tal]|2009.) . 

Data on ab undance variations w ithin the clusters are still rel- 
atively limited. iLetarte et al.l (1200 6*) noted that two of their nine 
stars show enhanced Na and depleted O abundances, as well as 
low Mg abundances. These patterns are similar to those observed 
in putative second-generation stars in Milky Way globular clus- 
ters. In our integrated-light measurements we observed Mg to 
be depleted relative to other a-elements (Ca, Ti) in Fornax 3, 4 
and 5, again suggesting that some stars in these clusters display 
anomalous abundances. It thus appears that the Fornax GCs are 
similar to the Milky Way GCs in this respect. 
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Fig. 1. Metallicity distributions for field stars and globular clus- 
ters in the Fornax dSph. The filled histog ram shows the raw dis - 
tribution of field star metallicities from iBattagha et al] (l2006h . 
while the outlined histogram is corrected for radial variations in 
the coverage (see main text). The hashed histogram represents 
the GCs with weights scaled by their luminosities and exagger- 
ated by a factor of 6 (note that Fornax 2 and Fornax 5 fall in the 
same bin). 



3. The globular cluster specific frequency revisited 

i Battaglia et al ] (I2006h measured metallicities for 562 field red 
giants in the Fornax dSph covering the full radial range out to 
~ 1°2, using spectroscopy of the Can near-infrared triplet ob- 
tained with the FLAMES spectrograph on the ESO VLT. About 
7% of the stars in their sample have metallicities [Fe/H] < -2. 
However, the raw sample may not accurately represent the true 
metallicity distribution of the Fornax stars as the metallicity dis- 
tribution does show a radial variation (with more metal-rich stars 
in the centre), and no correction has been made for the varying 
degree of completeness as a function of galactocentric distance. 
We have attempted to make a rough correction for this effect, 
proceeding as follows: 

We divided the profile of the Fornax dSph into 10 radial 
bins. For each bin, we compu ted the total number o f stars Nj(K) 
from the Sersic model fit of Battaglia et al. (2006) and counted 
the number of stars with FLAMES metallicity measurements, 
Nf{K). We then assigned a weight, W{K) = Nt{R)INf{R), to 
each star in the corresponding bin. Finally, a corrected metal- 
licity distribution was derived by counting each star accord- 
ing to its weight. The raw and corrected metallicity distribu- 
tions are shown in Fig. [T] Perhaps somewhat surprisingly, this 
procedure leads to a somewhat less pronounced low-metallicity 
tail, with about 5% of the stars having [Fe/H] < -2. This is 
because a much smaller fraction of the stars within a given 
FLAMES field are actually sampled in the central regions of 
the galaxy, so that the weights W(/?) increase towards the cen- 
tre where there are fewer very metal-poor stars. For comparison 
we also include the GCs in Fig.[T] Each cluster was "counted" as 
6 X A^stais X lo-o^CWi'.Gc+B.a) .^jjgj-g ^ is the total 
number of stars with FLAMES data, Mycc are the My magni- 
tudes of the individual clusters and My - -13.2 is the magni- 
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Fig. 2. Number of RGB stars (with My < - 0.5) per 10^ Mq 
formed versus age, based on Dotter et al.l (12007 ) isochrones and 
assuming a iSalpeteii d 19551) IMF normalized over the range 
0.15-100 Mo. 



tude of the Fornax dSph. The difference between the metallicity 
distributions of the field stars and the GCs is rather striking. 

If we simply scale the luminosity of the Fornax dSph by 5%, 
we obtain a total absolute magnitude of My = -10.0 for the 
low-metallicity tail. Relative to this tail, the four metal-poor GCs 
then correspond to a specific frequency of Sn = 400. With an 
absolute magnitude of My = -8.9 these GCs account for over 
1/4 of the integrated luminosity of all the stars with [Fe/H] < -2 
in Fornax! Note that this calculation assumes that numbers of 
RGB stars and integrated luminosities are equivalent, which is 
only approximately true as the number of RGB stars per unit 
luminosity depends somewhat on age and metallicity. 

Alternatively, we can sca le by stellar mass. F r om re solved 
colour-magnitude diagrams, IColeman & de Jon^ (12008 ) esti- 
mated that a total of 6.1 X lO^M© of stars have formed in 
Fornax (for a Salpeter IMF extending down to 0.15 M©). A 
so mewhat smaller mass of 4.3 x 10^ M© within 0?8 was found 
by Ide Boer et alj (l2012h . Although the star formation history is 
complex and extends up to a few 100 Myr ago, the vast majority 
of the stars formed more than 2-3 Gyr ago, and would therefore 
contribute to the RGB star population on which the Battaglia et 
al. metallicities are based. Converting the observed numbers of 
RGB stars per metallicity bin to stellar mass per bin is, how- 
ever, slightly non-trivial, as this conversion depends on both age 
and metallicity. Of interest here is the nu mber of RGB stars tha t 
are sufficiently bright to be included in the lBattaglia et al ] d2006h 
sample. The faintest stars in their sample are slightly fainter than 
V = 20, corresp onding to a limit of My =s -0.5. From stel- 
lar model tracks ( Dotter et al.ll2007l) . we find that the amount of 
time an RGB star spends at My < -0.5 is very nearly inde- 
pendent of age, but decreases with increasing metallicity from 
about 25 Myr at [Fe/H] = -2 to about 20 Myr at [Fe/H] = -1. 
If there are no age differences, low-metallicity stars would then 
be somewhat over-represented in Fig. [T] relative to the actual 
amount of mass present at low metallicities. To find the relation 
between number of RGB stars per unit stellar mass formed and 
age, the RGB lifetime must be multiplied by the specific evolu- 



tionary flux, essentially the product of the IMF (evaluated at the 
turn-off mass) and the rate at which the turn-off mass decreases 
(Renzini & Buzzoni 1986). The net result is that the number of 
RGB stars per unit mass decreases by roughly a factor of two 
over the time interval 3-12 Gyr at fixed metallicity (Fig. |2]i. 
This would then cause older stars to be underrepresented in 
Fig-El counteracting the trend with metallicity if there is an age- 
metallicity relation (Battaglia et al. 2006; Coleman & de Jon^ 
l2008h . If the bulk of the metal-rich stellar population is as young 
as 3^ Gyr and the metal-poor population is very old (> 10 Gyr), 
then Fig. |2] suggests that we will underestimate the mass of the 
metal-poor population by up to ~ 60%. The actual star forma- 
tion history has probably been somewhat more continuous over 
the range ~3-10 G yr (Battaglia et al. 2006; IColeman & de JongI 
120081; ldeBoeretah i2012), in which case Fig. [T| would more 
closely represent the true metallicity distribution by mass. From 
similar considerations, we find that the number of RGB stars per 
unit integrated luminosity varies much less with age, while the 
metallicity dependence is similar to that in Fig.|2l We may then 
have overestimated the fraction of light coming from metal-poor 
stars by -10-20%^ 

Scaling the IColeman & de Jon3 ("2008) mass estimate by 
5%, we then find that 3. 1 X IO'^Mq formed in stars with [Fe/H] < 
-2. About 35% of the stellar mass will h ave been lost due t o 
stellar evolution over a Hubble time (e.g. iLamers et al.ll2005h . 
leaving ~ 2 x 1O''M0 in metal-poor stars today. The combined 
masses of the four metal -poor GCs are ~ 1 -0 x lO^Mp assumin g 
an average M/Ly = 3.5 dPubath et al.lll992l;lLarsen et alJl2012h . 
From this we find that these clusters have a combined mass of 
about 1/2 of the mass in metal-poor field stars, and that about 1 in 
3 stars with [Fe/H] < -2 belongs to one of the four metal -poor 
GCs. If we correct for age effects, as discussed above, this ratio 
might decrease to 1:4 or 1:5. These estimates are very similar to 
that based on the luminosities. 

These are evidently rather rough estimates. The total mass 
of Fornax is uncertain, and the assumption that the metal- 
poor stars are smoothly distributed may not be valid if they 
formed i n a few discrete events associated with the GCs 
(Penarru b"ia"et al.l l2009b). Note that the FLAMES data prefer- 
entially sample fields close to, or including the GCs. Also, the 
dynamically based M/Ly ratio used for the GCs is higher than 
the M/Ly X 2 predicte d by SSP models for a standard IMF 
(e.g. IStraderetalJbOllh . Further, the Ca ii triplet-based metal- 
licity measurements may not be on exactly the same scale as our 
high-dispersion spectroscopy of the GCs, although this mainly 
affects the shape of the m etallicity distribution for [Fe/H] < -2 
dStarkenburg et al.ll2010h . Nevertheless, it is clear that a stag- 
geringly large fraction of the most metal-poor stars in Fornax 
belong to the four metal-poor GCs. For comparison, the Milky 
Way GC system has a combined mass of ~ 2.8 x lO^M© 
(iKruiissen & Porteg ies Zwartl l2009h . or about 2% of the mass of 
the stell ar halo if a t otal mass of ~ IO'Mq is assumed for the lat- 
ter (e.g. lSuntzeff et al.lll991D and we take into account that only 
~ 2/3 of the Galactic GC system is associated with the halo. 

4. GC formation and self-enrichment 

The fact that, after a Hubble time of dynamical evolution, a large 
fraction of the metal-poor stars in the Fornax dSph still belong 
to the GCs, has important implications for our understanding of 
the formation and dynamical evolution of GCs. 

First, there is little room for dissolution of any additional 
clusters formed together with the four extant metal-poor ob- 
jects. Each of these four surviving metal-poor GCs has likely 
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lost several times 10^ M p due to evaporation over a Hubble time 
(e.g. lJordan et al.ll2007h . so that these combined losses may well 
amount to over a million Mq. The implication is then that per- 
haps half of the most metal-poor stars now present in Fornax 
may originally have formed in these four GCs. This might sug- 
gest that the Fornax GCs did not form as part of a larger popula- 
tion of star clusters with a power-law mass distribution extending 
down to low masses. This is reminiscent of the star cluster pop- 
ulations in some present-day star forming dwarf galaxies which 
host a single, or a few, outstandingly massive clusters, one of 
the most extreme examples being NGC 1705 (iO'Connell et alj 
[19941 iHo & Filippenkol fl996). 

Second, these numbers represent a challenge for self- 
enrichment scenarios in which GCs were initially much 
more massive than they are no w Jp'Ercole et all lloost 
ISchaerer & Charhonne]ll201 ll:lBekkill2Ql Ih . It seems to be ruled 
out that the Fornax GCs could possibly have been more than 
a factor of a few more massive initially than they are now. 
This Umit may be even more stringent for Fornax 3, the most 
massive of the clusters, whose very low metallicity places it 
even further out in the tail of the field star metallicity distri- 
bution. Nevertheless, detailed abundance measurements suggest 
that the Fornax GCs share the same abundance anomalies seen 
in Galactic GCs. It has been suggested that the second gener- 
ation may have formed partly from gas that was subsequently 
accreted (Conrov & Spergel 201 1; Pflamm-Altenburg & Kroupi 
I2OO9I) . However, even this possibility does not easily explain the 
high ratio of GC to field stars in the Fornax dSph. It would re- 
quire that each cluster was able to accrete a significant amount of 
gas with the same composition as the first-generation stars, while 
little gas of this composition was available for the formation of 
additional field stars. 

An important assumption in the above discussion is that the 
field star population observed in Fornax today has not been sig- 
nificantly stripped by the Galactic tidal field. While this has 
likely affected some dSphs, the stellar component in Fornax it- 
self does not appear to have been affected by tidal stripping 
dPefiarrubia et alj|2009a ^. Conversely, it is possible that the GCs 
did not originate within the Fornax dSph, but formed in halos 
that were stripped of stars and dark matter before being accreted 
onto Fornax. This has also been suggested as a possible solution 
to the "timing problem", i.e., the fact that the clusters have not 
yet been dragged into the centre of the galaxy as expected due to 
dynamical friction (Cole et al. 2012). This scenario is difficult to 
rule out, but it is unclear why four or five globular clusters would 
have formed close enough to the Fornax dSph to be captured by 
it, given that no significant population of inter-galactic GCs has 
been found elsewhere in the Local Group. 

5. Conclusions 

We have used the most recent measurements of metallicities for 
globular clusters and field stars in the Fornax dwarf spheroidal 
galaxy to estimate the specific frequency of metal-poor GCs. 
Four of the GCs have [Fe/H] < -2.0; associating these GCs with 
field stars in the same metallicity range we estimate Sn ~ 400, 
by a wide margin the highest known anywhere. We find that 
about one quarter of all stars in the Fornax dSph with [Fe/H] < 
-2 belong to these four most metal-poor clusters. Accounting 
for standard dynamical evolution, an even larger fraction of the 
metal-poor stars must initially have been born in the clusters. 
This puts tight constraints on any further amount of stars that 
could have been lost from metal-poor star clusters due to early 
"infant weight loss" or "infant mortality". The extant clusters 



could at most have been a factor of ~ 4-5 more massive initially 
then they are now, and this requires a rather extreme scenario in 
which no field stars or other clusters of similar metallicity were 
formed initially. 

It would be of significant interest to carry out detailed chem- 
ical tagging of the most metal-poor stars in Fornax to determine 
how many of them can be traced back to the GCs. With sufficient 
accuracy, it should be possible to detect peaks in the field star 
metallicity distribution corresponding to stars lost from the GCs. 
One might also expect to see a significant fraction of metal-poor 
field stars with peculiar, GC-like abundance patterns. Fornax and 
other dSph galaxies may be unique laboraties in which the pres- 
ence or absence of anomalous abundances of light elements in 
the field stars might shed light on the origin of multiple stellar 
populations in globular clusters. 
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